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The Casimir-Polder (CP) potential and transition rates of a Rydberg atom above a plane metal 
surface at finite temperature are discussed. As an example, the CP potential and transition rates of 
a rubidium atom above a copper surface at 300k are computed. Close to the surface we show that 
the quadrupole correction to the force is significant and increases with increasing principal quantum 
number n. For both the CP potential and decay rates one finds that the dominant contribution 
comes from the longest wavelength transition and the potential is independent of temperature. We 
provide explicit scaling laws for potential and decay rates as functions of atom-surface distance and 
principal quantum number of the initial Rydberg state. 
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Rydberg atoms - atoms excited to large principal 
quantum numbers n - have attracted much attention in 
recent decades [l|. Aside from the inherent interest of 
studying such extreme states, the exaggerated proper- 
ties of these highly excited atoms make them ideal for 
examining the properties of a variety of systems that 
would be awkward to probe by other means. The large 
cross-sections and weakly bound outer electrons associ- 
ated with Rydberg atoms makes them extremely sensitive 
to small-scale perturbations and dispersion potentials, 
such as the van der Waals (vdW) and Casimir-Polder 
(CP) potentials 0. 

For example, the strong scaling of the free-space vdW 
potential between two Rydberg atoms with n (ex v}^) 
leads to the Rydberg blockade mechanism which has been 
put forward as a candidate for implementing controlled 
gate operations between isolated atoms [3|, Ij] . The effect 
relies on the massive level shift that one Rydberg atom 
experiences in close proximity to another. 

Level shifts of similar origin arise if the atoms are 
brought into the vicinity of a macroscopic body. With 
the increasing ability to trap and manipulate atoms close 
to macroscopic bodies, the effects of these surface (CP) 
potentials have become a subject of great interest. A 
plications ranging from novel atom trapping methods 
to atom chip physics [Gj. Thus, it is of both fundamen- 
tal and practical interest to understand the interplay be- 
tween atoms in highly excited states and field fluctuations 
emanating from macroscopic bodies. 

In this Letter we provide evidence that dispersion 
forces have a sizeable effect on the energy levels of highly 
excited Rydberg atoms when brought close to metallic 
surfaces, with shifts on the order of several GHz expected 
at micrometer distances. Due to the large atom size, 
next-to-leading order terms in the multipole expansion 
of the radiation field give additional contributions in the 
MHz range. Despite the existence of large numbers of 



thermal photons at 300K at the relevant atomic transi- 
tion frequencies, the level shifts are in fact temperature- 
independent [7[. 

For a given atom-field coupling iJjnt , the CP potential 
for an atom in state \n) and the radiation field in state 
\q) is given by the position-dependent part of the energy 
shift which, to second order in perturbation theory, reads 
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where En+q are the unperturbed energy eigenvalues of 
the atom-field system. In the long-wavelength approx- 
imation, the electric field couples to the atomic dipole 
moment d via the interaction Hamiltonian 



iJi„t = -d-E(r^), 



(2) 



with the electric field given in terms of the classical Green 
tensor (for a recent review see, e.g. [Si]) 
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The Green tensor G(r,r',a;) solves the Helmholtz equa- 
tion for a point source and contains all the information 
about the geometry of the system. The bosonic vec- 
tor fields fA(r, oj) describe collective excitations of the 
electromagnetic field and the linearly absorbing dielec- 
tric matter. 



The CP potential at temperature T acting on an atom 
in state \n) via a dipole interaction ([2]) is given by [9| 
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where • denotes the Frobenius inner product {A • B = 
^^ ^ y4.ij....ij,_B,;j^...ij.) and the primed summation means 
that the term with j — Q contributes only with half 
weight. Here G*^^^(r^, r^, w) is the scattering part of the 
Green tensor. The atomic polarizability is defined as 
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with ujkn = {Ek — En)/h denoting the atomic transition 
frequencies. The frequencies ^j = 2TTkBTj /h, i € N are 
the Matsubara frequencies and n{uj) = ^c^'^/'^bT _ j^j-i 
is the thermal photon number distribution. 

The reflective part of the scattering Green tensor of an 
infinitely extended planar metal that fills the lower half 
space z < is given by ^10|] 
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with r|| = (x, j/,0), k|| — {kx,ky,0) and k\\ = |k|||. The 
reflection tensor i?(k|| , z, z' , w) has the form 
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Here the unit vectors for s-polarized and p-polarized 



waves are e^ = ^ku x ^z and e^ 



p - (fc||e^T/3+efc||)/<7. 



The functions r^ = [£(aj)/3+ — /?_]/[£(aj)/3+ + /3_] and 
Tp — [/3+ — /3_]/[/3+ + /3_] are the usual Fresnel re- 
flection coefficients f or those wav es with wave numbers 



fcn, and q = Lo/c. 



P- = ^q^e[uj) - /c^ ana p+ = ^ q- 

The permittivity of the metal surface is modelled by the 
Drude relation, s{uj) = l — u!'^/uj{uj + ij), where Wp and 7 
are the plasma frequency and the relaxation rate of the 
metal, respectively. Magnetic effects will be neglected. 

Matrix elements of the dipole operator d = ef = efe^ 
for the transition between two electronic states |n, /, j, m) 
(n, principal quantum number; /, j, m, quantum num- 
bers for orbital and total angular momentum and z- 
component of the latter) and \n' ,l',j',m') factor into a 
radial and an angular part according to 

{n',l\j',m'\d\n,l,j,m) 

= e{Rn:^i,^,,\f\Rn,u){l'j'm'\er\ljm), (10) 



where \Rn,i,j) are the radial wavefunctions. The radial 
matrix elements are computed numerically using the Nu- 
merov method [ll|, ll2] in which the suitably scaled radial 
Schrodinger equation is integrated inwards until an in- 
ner cut-off point (commonly the radius of the rump ion) . 
The eigenenergies are computed as Enj.j = —Ry/n*^ 
(Ry: Rydberg constant) where n* = n — Sn^ij is the ef- 
fective quantum number and Sn^ij the quantum defect 
[l3| whose values are tabulated in the literature [lj|. 

To evaluate the angular part, we first convert from 
the jm basis to a mims basis {mi, tti^: ^-components of 
orbital angular momentum and spin) by summing over 
the relevant Clebsch-Gordan coefficients, 

(I'j'm'leJljm) = V Cii^'^/^ C''f-'V^ 
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with the orbital-angular momentum eigenstates lY/.m,) 
being spherical harmonics. Matrix elements in the mims 
basis are computed by rewriting the radial unit vector in 
terms of spherical harmonics [Y/m = yim{^, v)] 






and using the integral relation [dCl = sin i? d-d dip] 
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that expresses the angular integral over three spherical 
harmonics in terms of Wigner 3j-symbols. 

When Rydberg atoms are held sufficiently close to a 
surface, its effective radius (r) ~ aon^ [oq: Bohr ra- 
dius] can be on the order of micrometres and therefore 
a significant fraction of the surface distance. The dipole 
approximation is then no longer appropriate. In other 
words, the atom cannot be viewed as a point-like particle, 
and its non-negligible size requires the inclusion of con- 
tributions from higher-order multipoles. This correction 
can be found via a similar method as described above, 
with the dipole interaction Hamiltonian replaced by the 
quadrupole interaction Hamiltonian J15l | 
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In close analogy to the dipole case the CP potential for 
a quadrupole interaction is found to be 
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with the quadrupole moment operator Q = e(f (S) f )/2 
and the atomic quadrupole polarizability defined as 
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The matrix elements for the quadrupole transitions can 
again be evaluated by factoring Q — {e/2)Per ® e^ and 
computing the matrix elements between the radial and 
angular parts of the wave functions separately. Evalua- 
tion of the radial integral is again performed numerically. 
The tensor product of unit vectors in spherical harmonic 

form reads e^ ® e^ = 



If A with 
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The angular matrix elements can then be evaluated using 
Eq. dni). 

As can be seen from Eqs. dH) and (fTSJ) . the CP poten- 
tial is comprised of a pair of sums, one over the Matsub- 
ara frequencies and one over all available atomic tran- 
sitions. It turns out, however, that due to the finite- 
temperature environment only a limited number of dipole 
and quadrupole transitions contribute significantly to the 
total level shift. This effect is depicted in Fig. [1] where 
we show the relative contributions of the dipole transi- 
tions 43s — )• np and quadrupole transitions 43s — >■ nd to 
the total level shift of the 43s state of ®^Rb. Note that 
the final n state for the dominant dipole and quadrupole 
transition is shifted. This is due to the differing quan- 
tum defects for the p and d states. Moreover, for each 
of these individual (long- wavelength) transitions the first 
term in the Matsubara sum (with j = 0) dominates at 
the micrometer atom-surface distances envisaged here, 
and all other terms can be safely neglected. Remark- 
ably, we observe that the CP potential is independent of 
temperature from T = to 300K and beyond. As was 
recently shown [7J , this is due to the dominance of con- 
tributions from transitions whose wavelengths far exceed 
atom-surface separations. 

Figure[2^ shows the total CP potential (and hence level 



shifts) C/cp = C^cp + ^cp' ^'-'^ various ns states (with 
n = 32,43,54) of ^^Rb near a copper surface at 300K. 
As we are not interested in a particular transition chan- 
nel, the weighted sum over all possible final states has 
been taken. One observes that for very small (yet exper- 
imentally achievable and indeed desirable) distances of 
less than 2^m the expected level shifts rapidly grow to 
GHz sizes. At these distances, we also observe significant 
deviations of the total shift from the dipole contribution 
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FIG. 1: Relative contributions from different transitions to 
the CP dipole and quadrupole level shift of the state 43s of 
«^Rb. 
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FIG. 2: (a) Casimir-Polder level shifts of the states 32s (dot- 
ted line), 43s (dashed line) and 54s (solid line) of *^Rb near a 
copper surface at 300K: total shift and quadrupole contribu- 
tion alone (inset), (b) Spontaneous decay rate near a copper 
surface at 300K for the initial states 32s (dotted line), 43s 
(dashed line) and 54s (solid line) of *^Rb. 



([6]) alone due to the increasingly important quadrupole 
shifts [Eq. ([T5|) ] which themselves can be as large as sev- 
eral MHz (inset in Fig. Hi). 

Related to the energy level shift is a line broadening 
effect, i.e. an increased rate of spontaneous decay due to 
strong non-radiative processes, as the atom approaches 
the surface [8|, |l5| . This strongly enhanced body- induced 
spontaneous decay partially counteracts the expected in- 
crease in lifetime as a function of the principal quantum 
number n in free space (Fq oc n~^) [l|. 

In Fig. [2l3 we show the total decay rates of the Ryd- 
berg states ns {n = 32,43,54) of ^''Rb as a function of 
atom-surface distance. The body-induced decay rates for 
electric dipole and quadrupole transitions are calculated 
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Note that, unlike the CP potential, the decay rates are 
always temperature-dependent. One observes a strong 
increase of the decay rates near the surface (z < lO/im) 
which becomes more pronounced for states with higher 
principal quantum number n. This translates into a rel- 
ative line broadening of more than three orders of mag- 
nitude that potentially limits trapping and manipulation 
times of high- lying states near surfaces. For larger dis- 
tances (z > 15/^m) the rates quickly approach their free- 
space values and show the expected suppression with in- 
creasing n. 

Finally, we will briefly consider how the CP potential 
and transition rates scale with atom-surface distance z 
and the principal quantum number n. For metal surfaces, 
the reflection coefhcients are nearly independent of w at 
infrared frequencies and below. In the low-temperature 
limit, when the thermal photon number is negligible, 
there is no tu dependence for either the CP potential or 
the transition rates. The dipole and quadrupole moments 
for the dominant ns — >■ (n — l)p and ns -4- (n — l)d tran- 
sitions scale as n^ and n*, respectively, for large n. In 
the non-retarded limit (valid for surface distances beyond 
even 100 /^m), the body-induced rates and the CP poten- 
tial scale as z~'^ and z~^ for the dipole and quadrupole 
contributions, respectively il5i] . Combining these results 
leads to a scaling behaviour of 
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for the dipole and quadrupole components of the CP po- 
tential Ucp and the decay rate Tnk- 

In the high-temperature limit, the scaling of the 
CP shifts remains the same due to the temperature- 
independence demonstrated in Ref. (7| , whereas the tran- 
sition rates become proportional to the mean photon 
number n{uj) w ksT /{Hlo). For the dominant dipole 
and quadrupole transitions (Fig. [T]), one finds lu cx n~'^. 



and the transition rates scale as F^'^ cx riJ j z^ and 
^"^k °^ n^^/z^, respectively. 

We have shown in this Letter that the interaction 
between highly excited atoms and macroscopic surfaces 
leads to energy level shifts that can be as large as several 
GHz. This implies that any scheme that relies on the ma- 
nipulation of (trapped) Rydberg atoms near surfaces has 
to account for this major adjustment. Moreover, some of 
the advantages of using highly excited Rydberg atoms, 
in particular their rapidly decreasing Einstein coefficients 
with increasing principal quantum number n, are coun- 
teracted by the atom-surface interactions. 
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